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Express Mail Label No. EV 396082112 US 

[0001] SIMPLIFIED BLOCK LINEAR EQUALIZER 

WITH BLOCK SPACE TIME TRANSMIT DIVERSITY 

[0002] Cross Reference to Related Applications 

[0003] This application is a continuation of U.S. Patent Application Serial No. 

10/034,793, filed on December 27, 2001, which in turn claims priority from Provisional 
Patent Application No. 60/263,915, filed on January 25, 2001, which are incorporated 
by reference as if fully set forth. 

[0004] BACKGROUND 

[0005] The present invention relates to communication systems imploring code 

division multiple access (CDMA) techniques. More particularly, the present invention 
relates to a transmission diversity scheme which can be applied to a CDMA 
communication. 

[0006] Spatial diversity has been proposed for support of very high data rate 

users within third generation wide band code division multiple access systems. Using 
multiple antennas, the systems achieves better gains and link quality, which results in 
increased system capacity. Classically, diversity has been exploited through the use of 
either beam steering or through diversity combining. 

[0007] More recently, it has been realized that coordinated use of diversity can be 

achieved through the use of space-time codes. Such systems can theoretically increase 
capacity by up to a factor equaling the number of transmit and receive antennas in the 
array. Space-time codes operate on a block of input symbols producing a matrix output 
over antennas and time. 

[0008] In the past, space-time transmit diversity systems have transmitted 

consecutive symbols simultaneously with their complex conjugates. This type of 
system, though, may result in symbol overlap at the receiving end. The amount of 
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overlap is dependent on the length of the impulse response of the propagation channel. 
In time division duplex (TDD) mode, this symbol overlap will have to be accounted for 
in the joint detection receiver. The joint detector will have to estimate the overlapping 
transmitted symbols and their conjugates, resulting in an increase in complexity of the 
joint detection. 

[0009] In order to alleviate this increase in joint detection complexity, systems 

have been created which transmit two similar but different data fields. The first data 
field, having a first portion, Di, and a second portion, D2, is transmitted by the first 
antenna. A second data field is produced by modifying the first data field. The 
negation of the conjugate of D2, -D2*, is the first portion of the second data field and the 
conjugate of Di, Di*, is the second portion. The second data field is simultaneously 
transmitted by the second antenna. 

[00 10] Although this diversity transmission scheme reduces receiver complexity, 

receivers for this scheme are still very complex. Such receivers utilize two joint 
detection devices. Each joint detection device recovers the data field transmitted from 
one of the antennas individually. Such an implementation deals with cross 
interference between the two transmitted data fields by dealing with each antenna's 
transmission separately. As a result, each joint detection device treats the other 
antenna's transmission as noise. The symbols recovered from each joint detection 
device are combined using a decoder to determine d x and d 2 . A block diagram of this 
system is illustrated in Figure 1. The receiver in such a system has a high complexity 
due to the use of two joint detectors resulting in higher receiver expense. 
[0011] Accordingly, there exists a need for alternate receiver implementations. 

[0012] SUMMARY 

[0013] The present invention is a method and system for receiving data 

transmitted using block space time transmit diversity (BSTTD) in a code division 
multiple access (CDMA) communication system. The system comprises a transmitter 
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for transmitting a first data field using a first antenna and a second data field using a 
second antenna and a receiver. The receiver includes an antenna for receiving the first 
and second transmitted data fields, and a BSTTD joint detector which determines 
symbols of the first and second transmitted data fields using a minimum mean square 
error block linear equalizer model and an approximated Cholesky decomposition of the 
model. 

[0014] BRIEF DESCRIPTION OF THE DRAWINGS 

[0015] Figure 1 is a block diagram of a prior art communication system 

employing space-time transmit diversity. 

[0016] Figure 2 is a block diagram of a receiver in accordance with the preferred , 

embodiment of the present invention. 

[0017] Figure 3 is an illustration of matrix structures for approximation of Block 

Space Time Transmit Diversity (BSTTD) in accordance with the preferred 
embodiment. 

[0018] Figure 4 is a flow diagram of the block space time transmit diversity joint 

detection method in accordance with the preferred embodiment. 

[0019] DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

[0020] Figure 2 is a block diagram of a receiver 10, preferably located at a user 

equipment (UE), in a CDMA communication system in accordance with the preferred 
embodiment of the present invention. Although it is preferable to have the receiver 
located at the UE, the receiver 10 may be located at the base station and operating on 
uplink communications. The receiver 10 comprises a BSTTD joint detection device 
(BSTTD JD) 12, a channel estimation device 13 and an antenna 16. The antenna 16 of 
the UE receives various RF signals including a first and second communication burst 
from a transmitter. 

[0021] The first and second communication bursts comprise the first and second 

data fields, respectively as described above. The first data field includes the first 

-3- 



I-2-0181.2US 



portion Dl and the second portion D2; the second data field includes the negative 
conjugate of D2, -D2* and the conjugate of Dl, Dl*. A typical communication burst 
has the two portions of the data fields separated by a midamble. The burst also has a 
guard period at the end of it to allow for different times of arrival between bursts. 
Each data field of one communication burst is encoded as the first data field, Dl, D2. 
Each data field of the other communication burst is encoded as the second data field, - 
D2*, Dl*. The respective data fields are spread and a midamble included to produce 
the first and second communication bursts, respectively. Each of the communication 
bursts are transmitted by a respective first and second antenna in a RF signal to the 
receiver 10. 

[0022] The received RF communication signal including the first and second 

communication bursts is demodulated and forwarded to the channel estimation device 
13 and BSTTD JD 12. The channel estimation device 13 processes the demodulated 
signal and forwards the channel information to the BSTTD JD 12. 
[0023] The BSTTD JD 12 receives the demodulated signal including the first and 
second communication bursts and the channel information from the channel estimation 
device 13. Using the channel information and the spreading codes of the transmitter, 
the BSTTD JD 12 estimates the data symbols of the first and second data fields of each 
communication burst, Dl, D2, -D2*, Dl* and combines Dl, D2, -D2*, Dl* to recover 
the original data field D. 

[0024] In accordance with the preferred embodiment of the present invention, the 

BSTTD JD 12 estimates the data symbols of each of the received data fields utilizing a 
simplified minimum mean square error block linear equalizer (MMSE-BLE) based 
detector. The BSTTD JD 12 of the present invention operates in accordance with the 
following. A and B are block banded versions of the propagation matrices of channel 1, 
associated with antenna 1, and channel 2, associated with antenna 2, respectively. 
They are rewritten as a 2x2 block matrix as follows. 

_A 2i A 22 _ 
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[0025] Accordingly, the received signal model for block space time transmit 
diversity is expressed as Equation 1. 



2 J L 5 22 



0 



-B. 



21 



°. ~ Bi : 

^21 ^22 J 



Equation 1 



[0026] Since the length of the blocks is much longer than the channel delay 

spread, the interference between adjacent blocks, A21 and B21, can be ignored, and the 
received signal model can be simplified to Equation 2: 

Equation 2 
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[0027] In order to estimate the data blocks a MMSE BLE algorithm for BSTTD 
may be used. Using whitening matched filtering, the data blocks can be represented 
by Equations 3 and 4 below. 



( B 22 H r 2 ) 



Equation 3 



d wmf 2 =\i H r 2 -[B u H r x ) , 



Equation 4 



[0028] 



The MMSE-BLE output is represented as Equation 5. 
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Equation 5 



E is shown in Equation 2. a 2 is the mean noise variance and I is the identity matrix. 
[0029] In the single antenna BLE, the major complexity for block STTD is due to 

the matrix inversion, which is preferably implemented with an approximate Cholesky 
decomposition. The block matrix representation of the correlation matrix for Cholesky 
decomposition is written as Equation 6. 
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D = E"E + a 2 I = 



D U D 2l" 

D 2X D 22 



Equation 6 



Dn, D22 and D21 are per Equations 7, 8 and 9, respectively. 

D u = A u " A u + (b 22 h B 22 ) + <t 2 I Equation 7 

D 22 = B n H B u + (An" An ) + a 2 1 Equation 8 

D 2l = (a^BbJ - B U H A U . Equation 9 

[0030] The lower triangular matrix for the Cholesky decomposition D = GG* 1 is 

written as Equation 10. 



G = 



G„ 0 
L G 21 G 22 



Equation 10 



Equations 11, 12 and 13 are relationships between G11, G21, G22, Dn, D21 and D22. 

G n G u H = D u Equation 11 

G 2l G n " = D 21 Equation 12 

G 22 G 22 " = D 22 - G 21 G 2 " Equation 13 

[0031] The estimated symbol sequence can be obtained by solving the following 

triangular systems per Equations 14, 15, 16 and 17. 

Gu^x = d^j, Equation 14 

G 22 ih 2 = d^* - G 2X m, Equation 15 

Gi2 H d mmse 2 = m 2 Equation 16 

Gnd^ = m, -G 2l "d mm J Equation 17 
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[0032] In a single antenna system, one Cholesky decomposition is required. The 

use of a diversity antenna increases the complexity of the decoding of the symbols by 
requiring two Cholesky decompositions (Equations 11 and 13) and one forward 
substitution (Equation 12). This increases the complexity of a BSTTD system over a 
single antenna system by more than two times. Moreover, the BSTTD decoder of this 
system does not cancel the interference of the first sub-block to the second sub-block, 
which results in more error in the detection. 

[0033] The following describes reducing the complexity further. From the 

structure of the transmission matrix, A22 and B22 can be represented by the block 
matrix forms with An and Bn as follows. 



A, 



0 A 3 



and 



B 22 = 



B, 



0 B, 



Equations 18, 19 and 20 are relationships between An, A22, Bu and B22. 



Aj 2 A 22 - A n A n + 



0 0 

0 A 2 H A 3J 



Equation 18 



B ii" B 22 - B i\ h B n + 



0 0 

0 B"B y 



Equation 19 



B 22 - A u B u + 



0 



0 



0 A"B y 



Equation 20 



[0034] Those skilled in the art will realize that A \\ , A \ , B£ , B£ and A \\ , are 
the block Toeplitz matrices, but A" A,, , B" ,B U and A"B U are not because of the lower 
right sub-blocks in the last terms of the Equations 18, 19 and 20. 
[0035] Equation 4 , by substituting Equation 18, becomes Equation 21. 
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Dn=A 22 H A 22 ^(B 22 H 'B 22 )\a 2 I- 

Equation 21 is block Hermitian. The solution of the Equation 7 can be approximated 
by the repeated version of Cholesky decomposition by ignoring the last term, i.e., 
Equation 22. 

G n G" = D n Equation 22 

Dn is per Equation 23. 

Ai = A ?2 A zi + ( B ?2 B 2i)* + ° 2 I Equation 23 

Equation 22 is the block Toeplitz matrix approximation. Its complexity is equivalent to 
the approximated decomposition in the single antenna case. Those skilled in the art 
will recognize that the above equations result in an approximation of Gii, reducing the 
complexity of the BSTTD JD 12. 

[0036] Further reduction in the complexity of the BSTTD JD 12 can be found in the 
approximation of G22. From Equations 11 and 12, Equation 13 becomes Equation 24. 

G 22 G 22 = D 22 - D 2X D x ; l D 2X H Equation 

24 

With the assumption that norm(D 22 )» norm[D 2l D n ~ l D 2l H ), Equation 24 becomes 
Equation 25. 

^22^22 ~ D 2 2 

Moreover, from the Equations 8, 19 and 22, Equation 26 results. 

D 22 =Dn- 



0 
0 



0 

A 3 "A 3 



Equation 21 



0 0 

0 b"b 3 



Equation 26 
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[0037] Similar to the approximation of Gii above, the above solution can be 

approximated to the repeated version of Cholesky's decomposition by ignoring the last 
term, which results in Equation 27. 

G 22 = G* u Equation 27 

[0038] By this approximation, G22 and, hence, D22 (Equations 8 and 13) do not 

need to be computed explicitly. Therefore, the complexity of Cholesky decomposition 
with BSTTD becomes the same as the single antenna system. 

[0039] The major complexity of BSTTD over single antenna is associated with matrix 
G21 the Equations 12, 15 and 17. The number of complex operations in Equations 15 
and 17 is the same as the nonzero elements of G21. The less nonzero elements, reduces 
the complexity of Equations 15 and 17. One approach to reduce complexity is to 
assume G 21 = 0. However, this approximation introduces an error into the solution, 
which is typically not desired. 

[0040] Therefore, another approach to reduce complexity is to approximate G 21 in 
accordance with the following. From Equations 9 and 12, Equation 28 results. 

G 2l G n H = D 2l Equation 28 

D 22 is per Equation 29. 

D 21 = (A? 2 B 22 )* - B^A^ Equation 29 

[0041] Equation 29 results in a block Toeplitz matrix. Its general solution, 

though, is too complex to be readily implemented due to it's multiple forward 
triangular system solutions. However, it can be simplified using the following 
properties: 

Property 1: The matrix D 21 is skew-symmetric block Toeplitz, i.e., D 2l = -D 2l T . 
The diagonal terms of D 21 are always zeros. 
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Property 2: All the entries of D 21 are zeros except the elements in the last 
column or in the last row of the sub-block matrix. (See Fig. 2 (a)) 

Property 3: The matrix G 21 has a block Toeplitz structure. 

Property 4: The matrix G 21 is lower block banded with its bandwidth equal to 

(L • K a - 1) . (See Fig.2 (b)). L is the number of the non-zero blocks at the first row 

or column block. It is equivalent to the length of intersymbol interference plus 
one, i.e., ij — Lasi + 1, where Li S i = ceil(W/SF), W is the channel length and ceil(x) 
denotes the smallest integer larger than x. Ka is the total number of active 
codes (physical channel), e.g., Ka = K + 1 with K DCH in BCH timeslot. 

[0042] The complexity will be dramatically reduced using the above properties 

and the approximation to the block banded matrix with the same sub-block structure 
as in property 2 for D 2l . This approximated structure is shown in Fig.3 (c). Fig.3 (d) 
shows the exact G 21 with different collar scale from Fig. 3 (b). The computation of G 21 
will be simplified by the above properties as well as the following approximations: 

Approximation 1: G 21 is upper and lower block banded matrix with its 
bandwidth (L-K a -l). 

Approximation 2: G 21 has the same structure as D 2l . 
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[0043] 



With approximation 1, the simplified G 21 can be represented by: 
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The block matrix representations of the correlation matrix D 21 and lower triangular 
matrix G 21 are written as Equations 30 and 31. 
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Equation 30 



D 2l = 



-d T n 
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Equation 31 
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d,, equals- : dij and fij per property 1 and approximation 2 have the following 
structure. 



d = 

ij 



0 
0 

0 

x 



The solution of 



0 
0 

0 

x 



and f u = 



FG» = D 2l 



0 
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d 
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x 



\Ka 
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KaKa 



is obtained by computing the first block and first row block per Equations 32 and 33. 

fmSn =-df n , n = 1,2,..., L Equation 32 



fuSi = fu8" , n = 2,...,L Equation 33 

1=1 

A - \a {i \ and D = \d H \ with the above matrix structure and the lower triangular 
matrix G = [g^]^ satisfies the matrix equation AG H =D. is the number of 

dedicated channels (DCH) and K a = K d +1 is the total number of physical channels in 

the broadcast channel (BCH) time slot. The first K d element at the last column vector 

is obtained by the division of the complex number to the real number as per Equation 
34. 

a nKa = -^ 5fi - > n = K d - 1 • Equation 34 

8 KaKa 

[0044] The last row vector of matrix A 13 obtained by one forward 

substitution of size K a , which is represented by Equation 35. 
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_ U KaKa 



Equation 35 



[0045] In addition, the right hand side of Equation 33 contains matrix 

multiplications. Each matrix multiplication can be considered as K d + (K d + 1) 2 
complex multipliers due to the zero elements. 

[0046] The BSTTD algorithm is simplified using the above approximation as 

follows: 

Equations 
(3), (4) 
(23), (29) 
(22), (32), (33) 



Operation 

• Matched filter: 

• Correlation computation: 

• Cholesky decomposition: 



• Forward substitution per Equations 36 and 37: 

G 11™1 = d wm /l 

G U™2 =d wmf 2 -\ G 2l™l) 

• Backward substitution per Equations 38 and 37: 



Equation 36 
Equation 37 



G \\ d mmse2 ~ ™2 



U U d mmse\ = m \ "^21 d mmse2 



Equation 38 
Equation 39 



[0047] The preferred embodiment is described in conjunction with the flow chart 

of Figure 3. The received signal is modelled by ignoring the interference between data 
blocks, such as per Equation 2 (Step 401). The received vector is whitening matched 
filtered, such as per Equations 3 and 4 (Step 402). A Cholesky factor of the form of 
Equation 10 is determined for a MMSE BLE solution (Step 403). A sub-matrix of G, 
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Gn, is then calculated by calculating a Cholesky factor of a sub-matrix of D, Dn (of 
Equation 7), as per Equation 22 (Step 404). Another approximation of a sub-matrix of 
G, G22, using the complex conjugate of G11, G11*, per Equation 26 is calculated (Step 
405). Another sub-matrix of G, G21 is approximated as being an upper and lower block 
banded matrix using Equations 31 and 32 (Step 406). The symbols of the two data 
fields, dmmse l and dmmse 2 are solved using forward and backward substitution per 
Equations 35, 36, 37 and 38 (Step 407). The original transmitted data is then 
determined by decoding dmmse l and dmmse 2 using decoder 15 (Step 408). 

[0048] While the present invention has been described in terms of the preferred 

embodiment, other variations which are within the scope of the invention as outlined 

in the claims below will be apparent to those skilled in the art. 

# # # 
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